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ABSTRACT 

We have measured the kinematic profile of the early-type (E/SO) lens galaxy in the system 0047—281 
(z=0.485) with the Echelle Spectrograph and Imager (ESI) on the W.M. Keck-II Telescope, as part of the 
^ Lenses Structure and Dynamics (LSD) Survey. The central velocity dispersion is a = 229 ± 15 km s -1 , 

and the dispersion profile is nearly flat to beyond one effective radius (i? e ). No significant streaming 
motion is found. Surface photometry of the lens galaxy is measured from Hubble Space Telescope 
images. From the offset from the local Fundamental Plane (FP), we measure an evolution of the effective 
mass-to-light ratio of A log M/ Lb = —0.37 ± 0.06 between z = and z = 0.485, consistent with the 
observed evolution of field E/S0 galaxies. (We assume = 1, f2 m =0.3 and Ha=0-7 throughout.) 
| Gravitational lens models provide a mass of Me = (4.06 ±0.20) x 10 11 M Q inside the Einstein radius 

of R E = (8.70 ± 0.07) kpc. This allows us to break the de generacy between velocity anisotropy and 
y— | ' density profile, typical of dynamical models for E/S0 galaxies. We find that constant M/L model, even 

with strongly tangential anisotropy of the stellar velocity ellipsoid, are excluded at >99.9% CL. The 
total mass distribution inside Re can be described by a single power-law density profile, pt oc r~ 7 , 
with an effective slope 7' = 1-901q 23 (68% CL; ±0.1 systematic error). Two-component models yield 
an upper limit (68% CL) of 7 < 1.55(1.12) on the power-law slope of the dark-matter density profile 
and a projected dark-matter mass fraction of 0. 41(0.54)1^5 (l^og) (68% CL) inside Re, for Osipkov- 
Merritt models with anisotropy radius = oo(R e ). The stellar M*/L values derived from the FP agrees 
' well with the maximum allowed value from the isotropic dynamical models (i.e. the "maximum-bulge 

solution"). 

The fact that both lens systems 0047-281 (z=0.485) and MG2016+112 (2=1.004) are well described 
inside their Einstein radii by a constant M*/L stellar mass distribution embedded in a nearly logarithmic 
potential - with an isotropic or a mildly radially anisotropic dispersion tensor - could indicate that 
E/S0 galaxies underwent little structural evolution at Z < 1 and have a close-to-isothermal total mass 
distribution in their inner regions. Whether this conclusion can be generalized, however, requires the 
• ■ . analysis of more systems. We briefly discuss our results in the context of E/S0 galaxy formation and 

' cold-dark-matter simulations. 

Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: evolution 
— galaxies: formation — galaxies: structure 



1. INTRODUCTION 

To understand the evolution and internal structure of 
luminous and dark matter in early-type galaxies (E/S0), 
we have started the Lenses Structure and Dynamics (LSD) 
Survey. From an observational point of view, the LSD Sur- 
vey aims at measuring - using the Keck-II Telescope - the 
internal (stellar) kinematics of a relatively large sample 
of E/S0 gravitational-lens galaxies in the redshift range 
z = - 1 (see Treu & Koopmans 2002, hereafter TK02, 
for a description of the survey and its scientific ratio- 
nale). The stellar kinematic profiles are combined with 
constraints from a gravitational-lensing analysis - in par- 
ticular the mass enclosed within the Einstein radius - in 
order to break degeneracies inherent to each method in- 
dividually. In this way, we are uniquely able to constrain 



the luminous and dark mass distribution and the velocity 
ellipsoid of the luminous (i.e. stellar) component. 

Several of the main questions about E/S0 galaxies that 
we aim to answer with the LSD Survey are: (i) What is 
the amount of dark matter within the inner few effective 
radii (R e )l (ii) Does the dark matter density profile agree 
with the universal profiles inferred from CDM simulations 
(Navarro, Frenk & White 1997, hereafter NFW; Moore et 
al. 1998; Ghigna et al. 2001)? (iii) Is there a universal to- 
tal (luminous plus dark-matter) density distribution that 
well describes the inner regions of E/S0 galaxies? If so, is 
it isothermal, as observed in the local Universe and often 
assumed in lensing analyses? (iv) How does the stellar 
mass-to-light ratio of E/S0 galaxies evolve with redshift, 
and does it agree with the evolution of the stellar popu- 
lations of field early-type galaxies as inferred from Fun- 



1 Based on observations collected at W. M. Keck Observatory, which is operated jointly by the California Institute of Technology and the 
University of California, and with the NASA/ESA Hubble Space Telescope, obtained at STScI, which is operated by AURA, under NASA 
contract NAS5-26555. 
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damental Plane (hereafter FP) measurements (Treu et al. 
1999, 2001a, 2002, hereafter T99, TOla, T02; Kochanek et 
al. 2000; van Dokkum et al. 2001)? (v) Docs the struc- 
ture of E/S0 galaxies evolve between z = and 1, or is 
the evolution of the FP purely an evolution of the stellar 
population? (vi) Is the stellar velocity dispersion tensor 
isotropic or anisotropic? 

First results from the LSD Survey were presented in 
two recent papers (Koopmans & Treu 2002, hereafter 
KT02; TK02), where we combined our measurement of 
the luminosity-weighted stellar velocity dispersion of the 
lens galaxy in MG2016+112 at z = 1.004 with the mass en- 
closed by the Einstein radius as determined from gravita- 
tional lens models (Koopmans et al. 2002). A robust con- 
straint was found on the slope of the total density profile of 
the lens galaxy, i.e. 7' = 2.0±0.1±0.1 (random/systematic 
errors) for p t oc r~ 7 . In addition, we were able to deter- 
mine the stellar mass-to-light ratio and constrain the slope 
of the dark-matter halo, leading to a relatively simple self- 
consistent picture of the lens galaxy: an old and metal-rich 
stellar component embedded in a logarithmic (i.e. isother- 
mal) potential observed at a look-back time of ~ 8 Gyrs 
- remarkably similar to many present-day E/S0 galaxies. 
Constant M/L models were ruled out at a very high confi- 
dence level and a significant mass fraction (~75%) of dark 
matter was shown to be present inside the Einstein ra- 
dius of about 13.7 kpc. Unfortunately, given the faintness 
of the galaxy at z=l, only a luminosity-weighted velocity 
dispersion could be obtained and only minimal constraints 
could be set on the anisotropy of stellar orbits. 

In this paper, we present Keck and HST observations, 
a gravitational-lens model, and a dynamical model of the 
lens galaxy in 0047-281 at z = 0.485 (Warren et al. 1996, 
1998, 1999). This galaxy is bright and sufficiently ex- 
tended that we were able to measure a spatially resolved 
velocity dispersion profile, thus setting unprecedented con- 
straints on the orbital structure of a galaxy at a look-back 
time of ~ 5 Gyrs. 

The paper is organized as follows. In Section 2, we de- 
scribe archival HST observations and spectroscopic obser- 
vations with ESI on the Keck-II telescope. In Section 3, we 
analyze these observations and determine the luminosity 
evolution of the lens galaxy with respect to the local Fun- 
damental Plane. In Section 4, we present a gravitational 
lens model from which we determine the mass enclosed 
by the Einstein radius. In Section 5, we present a model 
for the luminous and dark-matter distributions of the lens 
galaxy. In Section 6, the results from the dynamical mod- 
els are presented. Section 7 summarizes and discusses the 
results. 

In the following, we assume that the Hubble constant, 
the matter density, and the cosmological constant are 
H = 65/i 6 5 kms^Mpc -1 , Ct m = 0.3, and fl A = 0.7, 
respectively. Throughout this paper, r is the radial coor- 
dinate in 3-D space, while R is the radial coordinate in 
2-D projected space. 



2. OBSERVATIONS 
2.1. Hubble Space Telescope Imaging 

2 Obtained as part of the CASTLeS Survey 



Wide Field and Planetary Camera 2 (WFPC2) images 
of the system are available from the HST archive, through 
filters F555W and F814W 2 . In particular, the system has 
been imaged for 9500s in F555W on the Wide-Field Cam- 
era, and 2700s in F814W on the Planetary Camera. 




Fig. 1— Left: HST F555W image of 0047-281. The images 
are designated A, B, C and D from the upper right going clockwise. 
Right: Same image with a smooth galaxy model subtracted. Note 
the almost continuous ring structure including images A, B, and C. 



The images were reduced using a series of iraf scripts 
based on the iraf package drizzle (Fruchter & Hook 
2002), to align the independent pointings and perform cos- 
mic ray rejection. A subsampled (pixel scale 0705) image 
was produced for the F555W image. The final "drizzled" 
image in F555W is shown in Fig.l. Note that for the red- 
shift of the source, z = 3.595, Ly-a is redshifted to 5589 A 
(Warren et al. 1998) and the F555W magnitude there- 
fore includes the bright Ly-a emission. The F814W image 
is significantly shallower than the F555W image and al- 
though adequate for measuring the structural parameters, 
it is not particularly useful for the lens modeling, since 
the multiple images are faint and their positions cannot 
be accurately determined. 

Surface photometry was performed on the F555W and 
F814W images as described in T99 and Treu et al. (2001b; 
hereafter TOlb). The galaxy brightness profiles are well 
represented by an i? 1 / 4 profile, which we fit - taking the 
HST point spread function into account - to obtain the ef- 
fective radius (R e ), the effective surface brightness (SB e ), 
and the total magnitude. The relevant observational quan- 
tities of galaxy G in 0047—281 and their errors are listed 
in Table 1. Note that errors on SB e and R e are tightly 
correlated and that the uncertainty on the combination 
logi? e — 0.325* B e that enters the Fundamental Plane (see 
Section 3) is very small (~ 0.015; see Kelson et al. 2000; 
TOlb; Bertin, Ciotti & del Principe 2002). In the right 
hand panel of Fig. 1 we show the HST image after removal 
of a smooth model for the lens galaxy. Notice the nearly 
circular structure of the lensed image configuration around 
the lens galaxy. 

Rest frame photometric quantities listed in Table 1 - 
computed as described in TOlb - are corrected for galactic 
extinction using E(B-V)=0.016 from Schlegel, Finkbeiner 
& Davis (1998). 

2.2. Keck Spectroscopy 

We observed 0047-281 using the Echelle Spectro- 
graph and Imager (ESI; Sheinis et al. 2002) on the 
W.M. Keck-II Telescope during four consecutive nights 
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(23-26 July, 2001), for a total integration time of 20,700s 
(7xl800s+3x2700s). The seeing was good (0"6 <FWHM< 
0'.'8) and three out of four nights were photometric. Be- 
tween each exposure, wc dithered along the slit to allow 
for a better removal of sky residuals in the red end of the 
spectrum. The slit (20" in length) was aligned with the 
major axis of the galaxy. The slit width of l'/25 yields 
an instrumental resolution of 30 km s _1 which is adequate 
for measuring the stellar velocity dispersion and removing 
narrow sky emission lines. The centering of the galaxy 
in the slit was constantly monitored by means of the ESI 
viewing camera (the galaxy was bright enough to be visible 
in a few seconds exposure) and we estimate the centering 
perpendicular to the slit to be accurate to < 071. 



Rcdshift (G) 
F8UW (mag) 
F555W (mag) 
$B e ,F8i4W (mag/arcsec 2 ) 
SB e ^ F555W (mag/arcsec 2 ) 
Re,F8UW (arcsec) 
Re Fbbbw (arcsec) 
b/a=(l - e) 
Major axis P. A. (°) 
a (kms -1 ) 

My - 5 log h 65 (mag) 
M B - 51og/i65 (mag) 
Re,v (h^kpc) 
R e ,B (h^kpc) 
SB e ,y (mag/arcsec 2 ) 
SB e ,s (mag/arcsec 2 ) 



0.485±0.001 
18.67±0.09 
20.61±0.09 
20.23±0.22 
22.59±0.07 
0.82±0.12 
0.99±0.04 
0.80±0.10 
67 ±5 



229±15 
-22.90±0.04 
-22.22±0.11 
5.21±0.72 
5.82±0.58 
19.29±0.22 
20.39±0.11 



Table 1. — Observed spectro- photometric quantities of the 
lens galaxy (G) in Q0047-281. The second part of the table 
lists rest-frame quantities, derived from the observed quanti- 
ties as described in Section 2. Note that a is the central veloc- 
ity dispersion corrected to a circular aperture of radius R e /8. 
All quantities in this table assume Ho = 65 kms _1 Mpc _1 , 
O m = 0.3, and D. A = 0.7. 

Data reduction was performed using the iraf package 
EASI2D 3 as described in KT02. To preserve most of the 
spatially resolved information and to achieve an adequate 
signal-to-noise ratio at the largest distance for the center 
of the galaxy, we defined 5 apertures along the spatial di- 
rection of the spectrum and summed the signal within each 
aperture. The apertures correspond approximately to an- 
gular dimensions along the slit: 076 x l'/25, Of! 6 x 1725, 
Of! 4 x 1725, 076 x 1725, 076 x 1725 (centered respectively 
at -171, -075,0", 075, 171 along the major axis). A ve- 
locity dispersion profile was measured with the Gauss- 
Hcrmitc Pixel-Fitting Software (van der Marel 1994) and 
the Gauss Hcrmite Fourier Fitting Software (van der Marel 
& Franx 1993) on the spectral region including the G- 
band (~ 4304 A), using as kinematic templates spectra 
of G-K giants observed at twilight with a 073 slit width, 
appropriately smoothed to match the instrumental reso- 
lution of the 1725 slit. The two codes provide consistent 
measurements. The total error on velocity dispersion was 
estimated by adding in quadrature the formal uncertainty 
given by the codes, the scatter in the results obtained with 
different templates and the semi-difference of the results 

3 developed by D. Sand and T. Treu; Sand et al. (2002), in prep. 



obtained with the two codes. The fit to the spectrum from 
the aperture including part of the brightest lensed images 
(A and B) was poor, with severe mismatch and unstable 
measurements. We interpret this as due to contamination 
by emission from the lensed images. We discarded this 
measurement from the analysis. The velocity dispersion 
profile was then folded around the center, as determined 
by fitting the centroid of the light distribution at the wave- 
length of the G band, and the velocity dispersion averaged 
in the corresponding apertures (symmetric apertures pro- 
vided results within the errors). The final results are listed 
in Table 2. 



Aperture (□") 


a (km s 1 ) 


Act (km s" 1 ) 


(0.0-0.2) x 1.25 


219 


12 


(0.2-0.8) x 1.25 


212 


14 


(0.8-1.4) x 1.25 


205 


13 



Table 2. — Kinematic data along the major axis of 0047—281. 
The adjacent rectangular apertures are indicated. The seeing 
was 077 during the observations. 

Using the procedure detailed in TOlb, the value in the 
central aperture can be corrected to standard central ve- 
locity dispersion of a = 229 ± 15 km s _1 within a circular 
aperture of radius R e /8 (see also Sec. 6.1). No evidence for 
significant streaming motions (e.g. rotation) was found, 
with an upper limit of 50 km s _1 relative radial velocity 
between the center and the outermost aperture. 

3. THE FUNDAMENTAL PLANE AND THE EVOLUTION OF 
STELLAR POPULATIONS 

Early-type galaxies in the local Universe occupy approx- 
imately a plane in the three-dimensional space defined by 
the parameters effective radius (logi? e ), effective surface 
brightness (SB e ) and central velocity dispersion (logcr), 

log Rq = a log a + (3 SB C + 7 FP (1) 

known as the Fundamental Plane (hereafter FP; Dressier 
et al. 1987; Djorgovski & Davis 1987). 

In recent years, it has been shown that a similar cor- 
relation between those observables exists in clusters out 
to z ~ 0.8 (e.g. van Dokkum & Franx 1996; Kelson et 
al. 1997; Bender et al. 1998; Pahre 1998; van Dokkum et 
al. 1998; J0rgensen et al. 1999; Ziegler et al. 2001). The 
observed evolution of the intercept 7fp of the FP with 
rcdshift is consistent with the expectations of passive evo- 
lution models for an old stellar population (rcdshift of for- 
mation Zf > 2; sec e.g. van Dokkum et al. 1998). No 
evidence for a dramatic evolution of the slopes a and (3 of 
the FP with redshift is found with the available data (see 
J0rgensen et al. 1999, Kelson et al. 2000, and TOla for 
discussion). The correlation is observed to be tight also in 
intermediate redshift field samples, although a faster evo- 
lution of the intercept is found in the highest redshift field 
samples (to z ~ 0.7) and interpreted as evidence for sec- 
ondary episodes of star formation in the field population 
at z < 1 (T02; see also T99, TOla, Kochanek et al. 2000, 
van Dokkum et al. 2001, Trager et al. 2000). 

Assuming that galaxy G in 0047-281 lies on a FP 
with slopes similar to those in the local Universe - and 
pure luminosity evolution wc can determine the offset of 
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its effective mass-to-iight ratio (M/L oc a 2 10^ OASB " / R e ) 
with respect to the local relation, which is related to the 
evolution of the intercept by A log M/L = — 0.4A7pp//3. 
As the local FP in the B band we adopt the relation 
found by Bender et al. (1998), i.e. a = 1.25, f3 = 0.32, 
and 7fp = —8.895 — log/150. In this way, we obtain 
A log Mj L b = -0.37 ± 0.06. The error is dominated by 
the observed FP parameters of galaxy G and dominates 
uncertainties on the local FP relation. In Fig. 2 we plot the 
evolution of the effective M/L for cluster and field E/SO 
galaxies as function of redshift (dashed and solid line; from 
T02), together with the value obtained for galaxy G (large 
filled square). The effective M/L evolution for galaxy 
G is consistent with what is observed for field galaxies, 
i.e. faster than for the cluster sample, possibly indicat- 
ing younger luminosity-weighted stellar populations (see, 
e. g., TOla, T02). As described and discussed in TK02, we 
can use this measurement to infer the stellar mass-to-light 
ratio (M*/Lb) of galaxy G assuming that 

log(M*/L s ) 2 - Iog(M./L B )o + A log(M/L B ), (2) 
where the second term on the right hand side of the equa- 
tion is measured from the evolution of the FP, and the 
first term on the right hand side of the equation can 
be measured for local E/SO galaxies. Note that Equa- 
tion 2 uses the non-trivial assumption that the stellar 
mass is a redshift-independent function of the combina- 
tion of observables used to define the effective mass a 2 R e 
(for a full discussion sec T02, and TK02). Using the lo- 
cal value of (7.3±2.1)/i 65 M q /Lb,© determined from data 
by Gerhard et al. (2001) as in TK02, we infer M*/L B = 
(3.1 ± 1.0)h 65 M /L B , for galaxy G. 




J I I I I I I I I I I I I I I I I I I I 

0.2 0.4 0.6 0.8 1 
Redshift 



Fig. 2. — Evolution of the M/L as inferred from the evolution of 
the FP. The solid and dashed lines indicate linear fits to the M/L 
evolution for field (T02) and cluster E/SO galaxies (van Dokkum et 
al. 2001) respectively. The large filled square at z = 1.004 indicates 
the M/L evolution of galaxy D in MG2016+112 (KT02) while the 
large filled square at z = 0.485 represents galaxy G in 0047—281. 

4. GRAVITATIONAL-LENS MODEL 

The Einstein radius (-Re) and mass enclosed by the Ein- 
stein radius, i.e. Me = M(< _Re) are quantities required 
in the dynamical models that will be discussed in Sec. 5. 
Both Re and Me are very insensitive to the assumed mass 
profile (see e.g. Kochanek 1991), especially for highly sym- 
metric cases such as 0047—281 (see Fig.l). For consistency 

- to ensure uniform modeling throughout the LSD survey 

- we model the lens galaxy as a Singular Isothermal Ellip- 
soid (SIE; Kormann et al. 1994; see also Chen, Kochanek 



& Hewitt 1995; Kochanek 1995; Grogin & Narayan 1996; 
Koopmans & Fassnacht 1999; Cohn et al. 2001; Munoz, 
Kochanek & Keeton 2001; Rusin et al. 2002). Although 
the assumption of a SIE mass distribution might not be 
accurate enough for some applications (see, e. g., Saha 
2000, Rusin & Tegmark 2001, Wucknitz 2002), we again 
stress that this choice of the mass profile does not signifi- 
cantly bias the determination of the quantities used in our 
analysis, i.e. Me and Re. The reason is that the image 
deflection angles in four-image systems like 0047—281 are 
not only nearly the same, but also that the deflection an- 
gle is only a function of the enclosed mass and therefore 
is little affected by either the radial mass profile or the 
ellipticity of the lens (see §3 in Kochanek 1991 for a more 
detailed discussion) . Constraints on the mass density pro- 
file obtained in Sec. 5 are therefore virtually independent 
of the choice of lens mass model. 



Image Ax (") Ay (") S(x, y) (mas) 



G 


=0.000 


=0.000 




A 


+ 1.048 


+0.726 


10,10 


B 


+0.896 


-0.802 


10,10 


C 


-0.126 


-1.165 


10,10 


D 


-1.011 


+0.263 


10,10 



Table 3. — Centroid image positions of 0047—281 with respect 
to the centroid of the primary lens galaxy (G). 



1 -1 




CM 

i 7 " 

- i i i i i - 

1 0-1-2 

Arcsec. 

Fig. 3. — A simple gravitational lens model of 0047—281. The 
cross-haired ellipse indicates the position, position angle and ellip- 
ticity of the lens galaxy. The thick solid contour is the critical 
curve, whereas the dashed contours indicate the time-delay surface. 
To guide the eye, the lensed images are indicated by two contours 
corresponding to circles in the source plane with 15 and 30-mas ra- 
dius. The subpanel shows a zoom-in of the caustics. The source 
position is indicated by a star. 

The centroids of the four lensed images are used as con- 
straints on the lens model (Table 3), assuming errors of 
10 mas, i.e. a fifth of a drizzled pixel size. We do not use 
the flux ratios in our analysis since they are not only dif- 
ficult to measure accurately for extended images, but also 
change as function of position along the arcs and could 
be affected by differential dust extinction (e.g. Surpi & 
Blandford 2002). The position and position angle of the 
lens-galaxy mass distribution are set equal to those of its 
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surface brightness distribution (Table 1) and the presence 
of an external shear is allowed for. Hence, there are six free 
parameters (i.e. the ellipticity and velocity dispersion of 
the lens-galaxy mass distribution, the source position, and 
the external shear), eight constraints (i.e. the four image 
positions) and consequently two degrees of freedom. The 
velocity dispersion of the SIE mass model is defined such 
that the mass enclosed by the critical curve is equal to that 
inside the Einstein radius of a Singular Isothermal Sphere 
with the same velocity dispersion (Kormann et al. 1994). 



i — i — i — i — i — i — i — i - 



0047-281 




0.5 1 

Radius (arcsec) 



1.5 



indication of a regular behavior in the physical properties 
of early- type galaxies, as discussed in Sections 5, 6 and 7 
(see also Kochanek 1994 and Kochanek et al. 2000). 

The adopted SIE velocity dispersion corresponds to a 
circular Einstein radius of Re = (8.70 ± 0.07)/i 65 1 kpc 
(i.e. 6» E = 1"34 ± O'.'Ol) and an enclosed mass of M E = 
4.06 x lO 11 /^ 1 M . The errors on Re and Me are cor- 
related (both depend on ctsie): f° r fixed Re one finds 
that SMe/Me = 2 (<5ctsie/csie)- This error corresponds 
to about 3% for 5a = 4 km s _1 , i.e. the difference be- 
tween the model with and without an additional compan- 
ion galaxy, which we adopt as systematic error (the ran- 
dom error on ctsie is a negligible 0.5%). By adopting a 
wide range of mass profiles and/or cllipticities the enclosed 
mass changes by < 4% for symmetric four-image systems 
like 0047-281 (e.g. Kochanek 1991). Hence, adding the 
two contributions in quadrature, the total error becomes 
5% on Me inside a radius of Re = 8.70h^ kpc. 

5. DYNAMICAL MODEL 

Following TK02, we model the galaxy mass distribution 
as a superposition of two spherical components, one for 
the luminous stellar matter and one for the dark-matter 
halo 4 . The luminous mass distribution is described by a 
Hernquist (1990) model 

I \ M * r * /Q\ 

PL(r) = 7r—r , - to , (3) 



27rr(r + r*) 3 ' 

where M* is the total stellar mass. This profile well 
reproduces the i? 1 / 4 surface brightness profile for r» = 
i? e /1.8153. In Sec. 6.1, we also examine the effect of a 
steeper inner core, Ph{f) oc r~ 2 , using the Jaffe (1983) 
model. We find the effect to be negligible within the er- 
rors. The dark-matter distribution is modeled as 



Pd(r) = 



Pdfi 



(r/r fc )7(l + (r/r h ) 2 )(3-7)/2 



(4) 



Fig. 4. — Velocity dispersion profile of 0047—281 along the major 
axis. The box height indicates the 68% measurement error, whereas 
the box width indicates the spectroscopic aperture (times the slit 
width 1'.'25). The open squares are the corresponding values for 
an isotropic constant M/L model, whereas the open stars are for 
a model with 10% lower Me (twice its error) and 10% higher R e 
(both of which lower a) with r-j = oo. The open circles indicate the 
isotropic Hernquist luminous component embedded in a total mass 
distribution with 7'=1.90. See Section 6 for details. 

The best solution for a single SIE mass distribution 
yields ctsie = 253 km s _1 for the lens galaxy, whereas the 
required external shear is 0.13. The x 2 is rather high, 306 
for two degrees of freedom. It is possible to device more 
elaborate models (for example with a dwarf companion) 
that fit the observational constraints much better. How- 
ever - as mentioned before - the mass enclosed by the 
Einstein radius is very insensitive to the precise lensing 
model. As an illustration of the uncertainty related to the 
mass modeling we considered including a dwarf companion 
galaxy with ctsie ~ 60 km s _1 . Although the fit is greatly 
improved (x 2 = because there are more free parameters 
than constraints), the velocity dispersion of the primary 
lens changes by only —4 km s~ 4 . We will therefore not 
elaborate further on the details of the lens model - which 
is not the intention of this paper - and use the results from 
SIE mass model in the rest of the analysis. More detailed 
lens models are being constructed based on the structure 
in the lensed arcs (R. Webster, private communication). 

We note that ctsie is close to the central stellar velocity 
dispersion (Table 1). However, we emphasize that ctsie 
is a model-dependent expression of the well-determined 
total enclosed mass, while the central stellar velocity dis- 
persion depends on the precise total mass profile, on the 
luminous mass profile, and on the dynamical state of the 
luminous component. Hence, the two quantities do not 
have to agree in principle. Their agreement is rather an 

4 Spherical dynamical models provide an adequate approximate description for computing kinematic quantities of E/S0 galaxies with axis ratio 
b/a ~ 0.8 like 0047-281 (c. g. Saglia, Bertin & Stiavelli 1992; Kronawittcr et al. 2000), even though clearly they are not appropriate for the 
lensing analysis, since, e. g., they do not produce quadruply lensed images. 



which closely describes a NFW profile for 7 = 1, and has 
the typical asymptotic behavior at large radii found from 
numerical simulations of dark matter halos oc r~ 3 (e.g. 
Ghigna et al. 2000). See TK02 for further discussion of 
this mass profile and dynamical model. 

According to the CDM simulations given in Bullock et 
al. (2001), a galaxy with the virial mass of galaxy G at 
z = 0.485 has w 50 kpc 3> R e ~ -Re- Hence, in the 
light of a comparison with CDM models, we can safely 
assume that the dark matter profile in the region of in- 
terest (i.e. inside ~ Re) is well described by a power-law 
Pd oc r~ 7 . Throughout this study, we will set r\> = 50 kpc 
(effectively equal to rt = 00). The effects of changing rt 
are discussed in TK02. 

In addition, we assume an Osipkov-Merritt (Osipkov 
1979; Merritt 1985a,b) parametrization of the anisotropy 
P of the luminous mass distribution 



r 2 _|_ r 2 ' 



(5) 



where ag and ay are the tangential and radial component 
of the velocity dispersion and r, is called the anisotropy 
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radius . Note that (3 > by definition, not allowing tan- 
gentially anisotropic models. A brief discussion of tangen- 
tially anisotropic models with negative constant values of 
(3 is given in Sec. 6. 2. As a further caveat, note that at 
infinite radii, Osipkov-Merrit models become completely 
radial. Although this behavior is not commonly found 
within the inner regions of E/SO probed by observations 
(e. g. Gerhard et al. 2001; but see van Albada 1982 and 
Bertin & Stiavelli 1993 for theoretical grounds), it has lit- 
tle impact in the case considered here, since the pressure 
tensor only becomes significantly radial well outside the 
Einstein Radius and in projection is significantly down- 
weighted by the rapidly falling luminosity-density profile. 

The line-of-sight velocity dispersion is obtained solving 
the three dimensional spherical Jeans equation (e.g. Bin- 
ney & Tremaine 1987 Eq. 4-30) for the luminous compo- 
nent in the total gravitational potential and computing the 
luminosity-weighted average along the line of sight (Bin- 
ney & Tremaine 1987; Eq. 4-60; see also, e. g. Ciotti, Lan- 
zoni & Rcnzini 1996). We correct for the average seeing 
of 0.7", during the observations, and average the velocity 
dispersion - weighted by the surface brightness - inside 
the appropriate rectangular apertures. For completeness, 
we rescale the apertures (Table 2) by 0.9 and 1.1 in the 
directions of the major and minor axes, respectively, such 
that their projection on the axisymmetric model is equiv- 
alent to their projection on an elliptical galaxy with an 
axial ratio of 0.8, even though this has minimal effects 
on the model velocity dispersions (< 1%), much smaller 
than the observational errors. The uncertainties on seeing, 
aperture size, and galaxy centering are taken into account 
as systematic errors in the following discussion. 

6. LUMINOUS AND DARK MATTER IN 0047— 281 

The unknown parameters of our dynamical model are 
M*, 7, Ti and pd,o (we note again that = i? e /1.8153 and 
rb = 50 kpc). We can eliminate one of these parameters 
using Me, the mass inside the Einstein radius, which is the 
most accurately known constraint on the lens mass model. 
We choose to eliminate pd.o- In addition, we transform 
M* into the stellar mass-to- light ratio M„/Lb, fixing the 
model luminosity exactly to the value 1.2 x 10 11 Lb.©. 
Hence, values of M*/ Lb that we derive from the dynami- 
cal model, bears an additional uncertainty of 11%, i.e. the 
observational error on Lb, whereas M*, which is used in 
the model calculation, does not have this error. We use 
the average R e = (5.52 ± 0.55)hg^ kpc of the rest-frame 
V- and B-band values. 

For any given set {M*/Lb,j} the dynamical model is 
completely determined and the luminosity- weighted veloc- 
ity dispersions for each of the three apertures (Tab. 2) 
can be computed. The likelihood is determined assuming 
Gaussian error distributions and the confidence contours 
using the likelihood ratio statistic 5 . 

6.1. Power- Law Models 

Before studying the luminous and dark-matter profiles 
individually, we determine the effective slope (7') of the 
total (luminous plus dark) matter density profile (p t ) in- 



side the Einstein radius. We emphasize that an effective 
slope of 7' does not imply that the density profile follows 
Pt oc r 1 exactly, but only effectively. 

For MG2016+H2 (TK02) we measured an effective 
slope 7' = 2.0 ± 0.1 ± 0.1 (random and systematic er- 
rors), assuming p t oc r~ 7 . Based on the data in Warren 
et al. (1998) and Kochanek et al. (2000), we found a similar 
effective slope for 0047—281. However, the errors were rel- 
atively large due to a lack of an extended kinematic profile 
and the large uncertainty on the velocity dispersion given 
by Warren et al. (1998). With the data presented here we 
can perform a more accurate measurement, and set con- 
straints on the anisotropy of the velocity ellipsoid. 

In Fig. 5(c), we show the likelihood contours of 7' ver- 
sus the anisotropy radius r i: based on our extended kine- 
matic profile of 0047—281. Three main conclusions can 
be drawn: (i) for a spherical isotropic stellar distribution 
function (n -> 00) 7' = 1.90 ± 0.05 (68% CL) (see also 
Fig.4); (ii) a lower limit of n/R e > 0.7 (68% CL) can 
be set, implying that the velocity distribution function is 
isotropic in the inner regions of the galaxy; (iii) indepen- 
dent of ri, one finds 7' = 1.901^ (68% CL). 

To assess systematic errors on 7', we varied Me and R e 
by the total uncertainties ±5% and ±10%, respectively. 
The value of 7', required to fit the data, changes by ±0.05 
and ±0.03, respectively. Other potential sources of errors, 
such as seeing, aperture corrections, aperture offsets, etc., 
were found to be negligible. We add both contributions 
and we conservatively round up to a total systematic er- 
ror of ±0.1, to account for all potential minor sources of 
error. 

Although the outer region of galaxy G is well fit by a 
Hcrnquist (1990) model, the inner region with pl(t) oc r _1 
is less constrained due to the finite resolution of the HST 
images. We therefore examined power-law models with 
steeper inner luminosity density profiles, pl(t) oc r~ 2 , 
using the model from Jaffe (1983). For both = 00 
and Ti = R e , we find that the stellar velocity disper- 
sions change by < 6 km s . In general, Jaffe models give 
slightly lower velocity dispersions for a fixed value of 7' 
and the best-fit models therefore yield slightly higher val- 
ues of 7' (few hundredths) compared with the Hcrnquist 
models. However, the differences are not significant, given 
the errors on dispersion profile, and therefore we conclude 
that our results are insensitive to the precise shape of the 
inner luminosity density profile. 

For completeness, we note that our model with 7' = 
1.90 results in a central stellar velocity dispersion of 
er=221 km s _1 inside R e /8, in excellent agreement with 
the empirically derived value of er=229±15 km s _1 in 
Sec. 2. 2, confirming the self-consistency of our models. 

In conclusion, the total mass distribution of galaxy G is 
well-matched by a single power-law density profile, that is 
isothermal (i.e. p t oc r~ 2 ) to within <~5%, and the velocity 
ellipsoid of the luminous component is isotropic at least 
inside ~ 70% of the effective radius. Note that this limit 
is generally consistent with the limits set by other phys- 
ical considerations. In fact, strongly radial orbits would 
result in radial instability (e.g. Merritt & Aguilar 1985 



5 Note that the likelihood ratio statistic is distributed as a \ 2 only asymptotically, hence the interpretation of the likelihood ratio contours 
as confidence contours is an approximation. However, because the confidence contours (i.e. Fig.5) on which we base our results only mildly 
deviate from true ellipses - i.e. the limiting case for large numbers of constraints this approximation should work well. 
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find that Osipkov-Merritt models with initial Jaffe density 
profile are unstable for n < 0.3r where r w i? e /0.763 is 
the Jaffe half-mass radius; for the effects of a dark matter 
halo see, e.g., Stiavelli & Sparke 1991) or negative values of 
the distribution function (e.g. Ciotti 1999, and references 
therein). 
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Fig. 5. — Panels (a) and (b) show the likelihood contours of the 
inner slope of the dark-matter halo (7) versus the stellar mass-to- 
light ratio (M„/Lb) for the lens galaxy in 0047—281, for n = 00 
and Ti = Re, respectively (the mass enclosed within the Einstein 
Radius is fixed to Mg, see Section 6 for details). The dashed lines 
indicate the likelihood contours from the dynamical model only, 
whereas the solid lines combine the constraints from the dynamical 
model and the FP. Panel (c) shows the likelihood contours of the 
effective slope (7') versus the anisotropy radius (rj). All models 
assume rf, = 50 kpc. See text for a more comprehensive description 
of the models. 



6.2. Constant M/L Models 

A stellar mass-to-light ratio of M*/Lb= (5.4 ± 
0.8)/i65 M q /Lb i0 is required to account for the mass Me, 
enclosed by the Einstein radius. This is larger than the 
value M*/Lb= (3.1 ± 1.0)/i65 M q /Lb iQ derived from the 
FP evolution. If we had no kinematic information this 
could only be interpreted as marginal evidence for dark 
matter inside the Einstein radius. 

However, the velocity dispersion profile changes the sit- 
uation dramatically: no model where mass follows light 
can be found to fit the data. For example, in Fig. 4 we 
show the velocity dispersion profile for an isotropic veloc- 
ity distribution and constant M/L model (open squares). 
The dispersion falls too sharply with radius and this model 
can be excluded at the > 99.9% CL. Even increasing R e 
by 10% and decreasing Me by 10% (twice its error) - both 
resulting in a smaller stellar velocity dispersion - the fit 
(open stars) can still be excluded at the >99.9% CL. Set- 
ting Ti = R e for the latter model (not shown) worsens the 
fit and the model can be excluded at > 99.9% CL. 

Strong tangentially anisotropic models for lens galaxies 
(see e.g. Romanowsky & Kochanck 1999), although prob- 
ably not very likely, can lead to flatter velocity disper- 
sion profiles, even if M/L is constant. We therefore tested 
models with constant negative values of (3. For (3 < —1.5, 
we find that the dispersion profile indeed becomes flat. 
However, the predicted stellar velocity dispersion is much 
higher (i.e. 240-260 km s _1 ) than the observed values 
and the model is excluded at > 99.9% CL. Models with 
(3 < —1.5 can only fit the data when Me is lowered by 
^30%. This is at least six times the error on Me and is 
incompatible with any acceptable lens model of 0047—281. 
We therefore conclude that mass does not follow light in 
the lens galaxy G of 0047-281, but that M/L must in- 



crease with radius. 

This is a key illustration that knowledge of the enclosed 
mass, Mg, inferred from the gravitational lens models, 
breaks the mass-anisotropy degeneracy. 

6.3. Stellar Mass and Dark-Matter Slope 

We now turn to the two-component models described in 
Section 5, to assess the density profile of the halo. In Fig. 5 
we show the likelihood contours (dashed lines) as a func- 
tion of the stellar mass- to-light ratio (M*/Ls) and inner 
slope of the dark-matter halo (7). We choose two repre- 
sentative values for the anisotropy radius, rj = R e and 
rj = 00. Note that both the mass within the Einstein Ra- 
dius (Me) and rj are fixed for all models in Figures 5(a) 
and 5(b). Hence, for each panel we effectively have two 
free parameters and three data points corresponding to 
the velocity dispersion profile (a fourth is given by the 
M*/Lb measured by the evolution of the FP, Section 3). 
We find that: (i) For M*/L s — > 0, the slope of the dark- 
matter halo 7 approaches ~2, consistent with the findings 
of Sec. 6.1 (panel c). (ii) For increasing values of M»/Lb, 
the total density profile and the velocity dispersion profile 
become steeper. Hence, 7 has to decrease to fit the data. 

When we include constraints on M*/Ls from the Fun- 
damental Plane (Sec. 3), the limits tighten on both M*/Lb 
and 7 (solid contours in Fig. 5). Remarkably, we find 
that the stellar mass-to-light ratio inferred from the FP 
agrees well with the maximum stellar mass-to-light ra- 
tio allowed by our isotropic dynamical models (as for 
MG2016+112'; TK02). Such a "maximum-bulge solu- 
tion" is in that sense equivalent to the "maximum-disk 
solution" for spiral galaxies (e.g. van Albada & Sancisi 
1986). We note however that the uncertainties are still 
considerable. In addition, the following limits are found: 
(i) M*/Ls=3.2+ 09 ^65 

M /L B;Q and 7 <1.55 (68% CL) 
for ri = 00 and (ii) M*/LB=2.5i. u 'g ftgs M /Le. q and 
7 < 1.12 (68% CL) for n = R e . Because the velocity dis- 
persion profile steepens with increasing radial anisotropy, 
a smaller 7 is required to retain a relatively flat dispersion 
profile. This explains the mild degeneracy between n and 
7- 

In the context of adiabatic contraction (AC) models 
(e.g. Blumenthal et al. 1986; Mo, Mao and White 1998), 
the initial slope (7$) of the dark-matter halo, i.e. before 
baryons assembled in the dark-matter potential well, is in 
general expected to be shallower than the observed DM 
slope (7). We find that the difference between 73 and 7 
is relatively small, because the stellar mass in the case of 
0047—280 is quite extended (i.e. large effective radius) and 
therefore affects the dark-matter slope less than found in 
MG2016+112 (TK02), for example. Even in the absence 
of AC, however, we find that % = 1.5 (Moore et al. 1998; 
Ghigna et al. 2001) is inconsistent with the results from 
0047-281 at the 90%(68%) CL for n = R e (oo), whereas 
the NFW profile (7* = 1; Navarro, Frenk & White 1997) 
is consistent at the 68% CL for rj > R e . However, any 
mechanism, including AC, that steepens the initial slope 
by more that A7 = 7 — 7$ «0.6, would imply that the 
results from 0047—281 are inconsistent with CDM simula- 
tions. 

7. SUMMARY & DISCUSSION 
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We have presented HST and Keck observations of the 
gravitational lens system 0047-281. In particular, HST im- 
ages have been used to measure the surface photometry of 
the lens galaxy G and to build a simple lens model of the 
quadruple-image system. Kcck-ESI data have been used 
to measure a spatially resolved velocity dispersion profile 
extended beyond the effective radius, with exquisite accu- 
racy (~ 5%). We have combined all these measurements 
to study the internal structure and dynamics of the lens 
galaxy at z = 0.485, finding the following: 

(i) The offset of galaxy G from the local Fundamental 
Plane, A log Mj L b = -0.37 ± 0.06 between z = and 
z = 0.485, is consistent with what is observed for field 
E/S0 galaxies at similar redshift (T02), i.e. somewhat 
larger than for cluster E/S0 galaxies. In terms of pure 
luminosity evolution this could be explained with inter- 
mediate age single stellar populations, or - more likely - 
with secondary episodes of star formation contributing a 
fraction of young stars to an old underlying stellar popu- 
lation (see discussion in T02). 

(ii) No dark-matter or constant M/L models are excluded 
at > 99.9% CL. Also constant M/L models with strongly 
tangential anisotropy of the stellar velocity ellipsoid are 
excluded at >99.9% CL. 

(iii) The stellar mass-to-light ratio M*/L B =(3.1 ± 
1.0) /i65 M q /Lb !Q obtained from the offset of the FP is 
inconsistent with the required Mh»/L b =5.4 M /L B!0 
to fully account for Me- This suggest the presence 
of dark matter. The FP value is consistent with 
what is obtained with our two-component dynamical 
models, and combining the two constraints we find 
M*/L B =3.2j:g;|/i 65 M /L B , (68% CL) for an isotropic 
velocity ellipsoid. Hence, dark matter comprises a fraction 
of 0.4llg;°g of the total mass enclosed by the Einstein ra- 
dius of 8.70 hg^ kpc for rt — oo. For r,- L = R e this fraction 
increases to 0.54+QQg. The data are also consistent with an 
Osipkov-Mcrritt (OM) radial anisotropy with anisotropy 
radius n > 0.7 R e (68 % CL). 

(iv) The total (luminous plus dark) mass distribution in- 
side the Einstein radius can be described by a single power- 
law density distribution, p t oc r~ 7 , with 7' = 1.90 ± 0.05 
(68% CL) for isotropic models, i.e. ri = 00. In general, 
i = 1.90±g;^l (68% CL) is found. The systematic error is 
estimated at 0.1. 

(v) An upper limit 7 < 1.55 (68 % CL) is found on the 
slope on the dark-matter halo inside the Einstein radius 
for an isotropic model. This limit tightens to 7 < 1.12 
for mildly anisotropic models with j-j = R e . Initial dark- 
matter profiles with 7$ = 1.5 (Moore et al. 1998; Ghigna 
et al. 2001) are therefore only marginally acceptable, es- 
pecially since the profile is expected to be less steep before 
the galaxy assembled. If 7 steepens by A7 > 0.6 dur- 
ing galaxy formation all CDM simulations are inconsistent 
with our results. 

In summary, the lens galaxy in 0047—281 appears to 
convey the same picture formulated for MG2016+112 that 
early-type galaxies at significant look-back times can be ef- 
fectively described by a i? 1 / 4 luminous component (mod- 
eled in this paper as either a Hernquist or Jaffe profile) em- 
bedded in a nearly-isothermal total mass distribution and 
that their stellar velocity dispersion is relatively isotropic, 



in particular inside the effective radius. In fact, both lens 
galaxies in MG2016+112 (KT02, TK02) and 0047-281 
show that the total mass distribution is well approxi- 
mated to within 5% by a simple power law density profile 
Pt oc r~ 2 (i.e. isothermal). Whether this conclusion can 
be generalized, however, requires the analysis of more sys- 
tems. 

Even more so, we have shown that deviations from 
isothermality or isotropy in the lens galaxies of 0047—281 
and MG2016+112 quickly lead to inconsistencies with con- 
straints from either the FP, the observed stellar kinemat- 
ics, the stellar mass-to-light ratio, observations of local 
E/S0 galaxies, the gravitational-lens models, etc., whereas 
the models that fit all constraints are internally consistent, 
appear to agree with all observational constraints avail- 
able, and indicate both isothermality and near-isotropy. 
Constant M/L or steep mass profiles inside the Einstein 
radius are excluded at very high confidence levels. 

A physical explanation is required if isotropy and the 
almost perfect isothermality are confirmed to be generic 
features of early-type galaxies. In particular, this regu- 
larity might suggest that luminous and dark-matter were 
strongly coupled at some point during galaxy assembly. 
Whereas adiabatic contraction has been suggested as a 
mechanism that can lead to near-isothermal mass profiles 
(e.g. Keeton 2001), it is not clear why such a process should 
only stop when the density profile is isothermal to better 
than apparently a few percent (see also TK02). Adiabatic 
contraction also leads to a slope of the inner density profile, 
inconsistent with the observed absence of lensed images in 
the centers of lens galaxies (e.g. Keeton 2001), if either the 
central black holes are not very massive or the inner den- 
sity profiles do not steepen through some other process. 
Violent relaxation could be a natural and viable explana- 
tion for this regularity, although it should also be explained 
why luminous and dark matter have different density pro- 
files (see discussion in TK02 and references therein). A 
combination of the two processes during some period in 
the formation of the galaxy can not be excluded. 

The striking similarity of the internal structure of E/S0 
galaxies at large look-back times with the internal struc- 
ture of local E/S0 galaxies (e.g. Franx et al. 1994; Bertin 
et al. 1994; Rix et al. 1994; Gerhard et al. 2001; see also 
Kochanek 1995 and reviews by de Zeeuw & Franx 1991; 
Bertin & Stiavelli 1993; Merritt 1999) suggests that lit- 
tle structural evolution occurred during the past 8 Gyrs 
(although again a larger sample is needed to make a quan- 
titative and general statement). The lack of significant 
structural evolution is also suggested by the remarkable 
agreement between the stellar M*/L B obtained with our 
dynamical models and the stellar M*/L B estimated us- 
ing the FP evolution. This fact adds further evidence in 
favor of a scenario where the general population of mas- 
sive (field) E/S0 galaxies changed little in the past 8 Gyrs 
(from z ~ 1) — as indicated for example by the modest 
evolution in their number density (Schade et al. 1999; Im 
et al. 2002; Cohen 2002; McCarthy et al. 2002) and by 
the little evolution in the scatter of the FP (T02) - with 
most of the evolution being driven by ageing of old stars 
and secondary episodes of star formation (Jimenez et al. 
1999; Trager et al. 2000; Menanteau, Abraham & Ellis 
2001; T02). 
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